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Summary 

Nomographic  design  procedures  are  developed  for  a  metallic  cylinder  (eg 
a  rocket  motor  case)  circumferentially  reinforced  with  a  pre-stralned 
fibre  overwind  with  viscoelastic  properties  .  Techniques  for  calculating 
the  pressure  for  first  yield  in  the  metallic  case  and  the  optimum  winding 
conditions  for  both  short-  and  long-  times  are  presented,  together  with  a 
technique  for  determining  the  maximum  fibre  stress  . 


Design  examples  are  presented  with  each  nomogram  in  turn  and  the  graphical 
results  are  compared  with  exact  numerical  results  . 
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1  INTRODUCTION 

Recent  work  (Refs  1-4)  on  the  stress  states  in  metallic  cylinders  circum¬ 
ferentially  reinforced  with  a  prestrained  fibre  overwind  has  involved  analyses 
of  the  initial  winding  process,  the  effects  of  different  storage  temperatures, 
pressurisation  and  viscoelastic  relaxation  of  the  fibre  overwind  .  The 
work  has  particular  relevance  to  the  rocket-motor  case,  and  much  attention 
has  been  given  to  the  effects  of  the  material  and  design  variables  on  the 
pressure  required  for  the  initiation  of  yield  in  the  reinforced  cylinder  of 
such  a  case  . 

The  range  of  choice  through  the  several  design  variables  is  considerable, 
and  it  is  not  easy  to  assess  directly  the  relative  Importance  of  the  effects 
of  changes  in  these  variables  .  It  is  in  these  circumstances  that  the  use 
of  nomograms  for  the  representation  and  interpretation  of  the  design 
equations  offers  very  considerable  advantages,  which  far  outweigh  the 
small  Inaccuracies  usually  associated  with  graphical  solution  techniques  . 

In  this  report  design  procedures  are  developed  for  the  stress  analysis  of 
fibre-reinforced  metallic  rocket-motor  cases  .  These  procedures  are  pre¬ 
sented  graphically  in  such  a  way  that  the  effects  of  changes  in  the 
variables,  and  their  relative  importance,  can  be  easily  seen/  It  should 
be  noted,  however,  that  the  procedures  presented  are  general,  in  that 
they  may  be  applied  to  any  fibre-reinforced  cylinder  subjected  to  internal 
pressure  .  The  effects  of  the  overwind  variables  (including  viscoelastic 
relaxation)  and  pressurisation  are  taken  into  account,  but  thermal  effects 
have  been  omitted  .  Experience  has  shown  that  the  small  changes  in  the 
yield  pressure  resulting  from  typical  temperature  changes  do  not  warrant 
the  extra  complexities  in  the  design  charts  . 

Two  of  the  five  nomograms  presented  are  concerned  with  the  determination 
of  the  pressure  for  the  Initiation  of  first  yield  in  the  metallic  cylinder, 
two  with  an  optimisation  analysis,  and  the  final  nomogram  evaluates  the 
maximum  fibre  stress.  In  Che  last  nomogram  viscoelastic  effects  have 
been  neglected,  since  numerical  results  have  shown  conclusively  that  the 
maximum  fibre  stress  occurs  immediately  after  winding  in  the  outermost 
layer  . 

Design  examples  are  presented  with  each  nomogram  .  The  accuracy  of  the 
solutions  are  assessed  by  comparing  graphical  results  with  exact  numerical 
results  . 

2  ELASTIC  ANALYSIS 

2  .1  Yield  Initiation  Predictions 

2  .1 .1  Theory 

It  is  a  requirement  that  the  rocket-motor  case  must  withstand  all 
operational  loads  without  failure  or  excessive  distortion.'  Although  aero- 
heat  and  inertia  effects  may  contribute,  the  predominant  load  will  always 
be  the  firing  pressure  .  The  maximum  permissible  pressure  is  defined  as 
that  for  the  onset  of  first  yield  in  the  metallic  case.  Although  such  a 
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definition  does  not  take  into  account  the  strength  reserves  available  in 
the  case  after  the  initiation  of  yield,  this  definition  will  be  adopted 
here  • 


It  has  been  shown  (Ref  3)  that  for  a  thin-walled  overwound  motor  case  of 
radius  R,  overwound  with  n  elastic  fibre  layers  with  an  initial  winding 
strain  and  subjected  to  the  applied  firing  pressure  p,  the  resulting 
circumferential  stresses  in  the  l'th  fibre  layer  and  the  case, 
and  Ogj.p'*''  respectively,  are  given  by 


aefp(1’n)  ■  c„Ef  a-AFi(A,n)}  -  {1-  £  (1-^!)}  • 

1*1 ,2 ,  •  e  •  jH 

and 

°0cp(n)  “  -Ec6wAFl(A,n)  +  (1+  ^  Q-^>>  , 

where 


and 


A 


Ef*f 

Ectc 


F^A.n) 


I  _ 

j-i  1+Aj 


(1) 

(2) 

(3) 

(4) 


In  the  above  equations  Ef  and  Ec  denote  respectively  the  Young's  modulus 
of  the  fibre  and  case,  tf  the  effective  thickness  of  a  single  fibre  layer 
and  tc  the  thickness  of  the  case.  Re  the  effective  nozzle  radius  (Ref  2.) 
and  the  Poisson’s  ratio  of  the  case  material-  The  corresponding 

axial  stress  component  in  the  case  follows  from  the  transverse  equilibrium 
equation,  le 


(R^V) 

<»zcp  "  P  -  2Rt c  ' 


(5) 


Once  the  stress  components  in  the  metallic  case  (Eqns  2  and  5)  are  known, 
the  pressure  for  the  onset  of  yield  can  be  calculated  using  the  von  Mlses 
yield  criterion,  the  most  appropriate  for  the  case  materials  under 
consideration  .  If  the  radial  stresses  are  assumed  negligible  relative  to 
the  circumferential  and  axial  stresses,  the  yield  criterion  (Ref  5) 
reduces  to  the  form 


<ff9cp<n>>2  +  <*~n>2  "  °ae0(n) 


zcp' 


8cp 


"zcp 


(6) 
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where  Oy  is  the  uniaxial  yield  stress  of  the  case  material  •  Equations 
2,  S  and  6  form  the  basis  for  nomogram  construction. 

If  used  in  their  present  dimensional  form,  these  equations  would  result 
in  nomograms  of  unacceptable  complexity.  However,  if  the  equations  are 
changed  into  a  non-dimensional  form  the  number  of  variables  Is  reduced 
and,  as  a  consequence,  their  graphical  representation  is  considerably 
simplified  . 

The  following  non-dimensional  variables  are  introduced: 


6  “  An 

(7) 

X  -  l-R^/R2 

(») 

P/p 


ymax 


p/3Rl 

4tc°y 


+w  “  Ecew/oy 

«0w  “  <W°y  "  *wAF1(A.n) 


(9) 

(10) 

(ID 


* 

where  0&,  is  the  circumferential  case  stress  due  to  winding, 

®9cp<n)  “  °8cp(n)/°y 
°zcp  “  °*cp/°y 


(12) 

(13) 


In  Equation  9  the  quantity  Pymax  denotes  the  maximum  attainable  yield 
pressure  for  the  reinforced  motor  (Ref  2),  as  given  by  the  equation 


Pymax 


(14) 


For  analytical  purposes  it  is  convenient  to  introduce  the  approximation 


AFi(A,n)  ♦  ln(l+An) 


(15) 


so  that  Equation  11  can  be  expressed  in  terms  of  the  variable  B,  that  is 
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®0w  -  *wln(l+8) 


(16) 


With  these  non-dimensional  variables,  Equations  2,  5  and  6  can  be  expressed 
in  non-dimensional  form  as 


A 


+ 


2p 

✓*X(1+B) 


{2+  I*) 

3 


zcp 


2 

73 


A 

P 


(»ecP(n))2+  “ 


zcp' 


36cp 


(n)S 


zcp 


(17) 

(18) 

(19) 


where,  in  Equation  17,  it  has  been  assumed  that  vc  •  1/3,  a  typical 
value  for  most  suitable  case  materials  • 

The  variables  B  and  X  are  now  combined  by  introducing  the  auxiliary 
variable  $,  where 


♦  -  (2+  3-}/{x(i+b:: 


(20) 


so  that  Equation  17  reduces  to  the  four-variable  equation 


-a 


8w 


(21) 


Once  the  governing  equations  have  been  obtained  in  non-dimensional  form 
(Eqns  18,  19^and  21),  in  order  to  predict  the  pressure  for  the  onset  of 
gleld  (ie  p“p„)  it  is  necessary  to  eliminate  the  quantities  and 

°zcp  from  Equation  19  in  favour  of  Py  using  Equations  21  and  Its 
respectively,  to  give  the  quadratic 


(22) 


Once  py  has  been  obtained,  the  yield  pressure  py  readily  follows  from 
Equation  9  . 


8 

UNLIMITED 


UNLIMITED 


Equation  22  can  now  be  represented  graphically  by  plotting  the  variation 
of  ogu  with  p  for  various  py  values,  typically  from  0  to  1  in  increments 
of  0.1.  This  is  depicted  in  the  central  portion  of  Figure  1.' 

This  graphical  technique  for  evaluating  py  is  in  terms  of  the  variable 
combinations  ogw  and  4.  It  is  now  necessary  to  extend  the  Equation 
22  plot  to  enable  the  evaluation  of  py  to  be  carried  out  In  terms  of  the 
more  fundamental  variables  py,  8  and  X.  From  Equation  16  it  can  be 
seen  that  the  variable  ogy  is  a  simple  multiplication  of  py  and 
ln(l+6);  this  equation  can  therefore  be  represented  by  the  'z'  type 
nomogram  (Ref  6)  shown  to  the  left  of  the  central  plot  in  Figure  1  . 

By  contrast,  the  graphical  determination  of  p  necessitates  a  rearrangement 
of  the  governing  Equation  20  so  that  it  complies  with  an  alignment  chart 
(Ref  6)  of  the  form 


f l(X)  +  f2(S)f3(p)  -  *4(8) 

where  the  functions  represent 
f3(X)  -  -2/X 
f2(B)  -  1  +  8 
£  3(  ♦)  “  ♦ 
f 4( 6)  -  8/3 


(23) 


(24) 


The  Equation  23  nomogram  is  shown  in  the  lower  portion  of  Figure  1  . 

The  range  of  the  variables  associated  with  these  nomograms  must  be 
chosen  with  care  to  ensure  that  the  final  design  chart  can  be  used  for 
all  practical  applications  .  For  fibre-reinforced  rochet  motor  cases 
these  ranges,  together  with  those  to  be  used  in  subsequent  nomograms,  are 
presented  in  Table  1  . 

Since  the  variable  ogy  is  common  to  both  Equations  16  and  22,  and, 
furthermore,  since  the  variable  p  is  common  to  Equations  20  and  22,  it 
follows  that  the  three  nomograms  developed  can  be  combined  in  such  a  way 
that  the  transfer  of  data  from  one  nomogram  to  another  is,  in  general, 
avoided.  This  reduces  reading  errors  and,  omitting  constructional  details, 
results  in  the  simple  design  chart  presented  in  Figure  1 . 

There  are  four  fundamental  variables  in  the  design  process  embodied  in 
Figure  1,  py,  X,  py  and  8>  Any  one  of  the  first  three  (le  (^or  X 
or  py)  can  be  readily  'read'  from  the  Figure  1  design  chart  provided  the 
other  three  are  known.  However,  because  the  variable  8  Is  represented  in 
two  of  the  graduated  loci  in  the  design  chart.  It  cannot  be  obtained 
directly  In  terms  of  the  variables  py,  X  and  Py.  Nevertheless,  a  trial- 
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and-error  graphical  approach  for  g  Is  feasible  .  The  problem  Is  resolved 
If  Ofy  (le  the  case  stress  Induced  by  winding)  can  be  specified  together 
with  py  and  X;  In  these  circumstances  the  solution  proceeds  as  follows: 

a  .  Determine^the  4  value  corresponding  to  the  Intersection 

of  the  specified  ogu  and  py  values  in  the  ogw  ~  4  grid  depicted 
in  the  central  portion  of  Figure  1 . 

b.  Determine  the  required  g  value  from  the  Equation  23 
nomogram  (lower  portion  of  Fig  1)  using  the  calculated  4  and 
specified  X  values  . 

c  .  Obtain  ^  and  thence  the  winding  strain  which  satisfies 

the  prescribed  problem  using  the  left-hand  side  of  Figure  1  . 

2  .1  .2  Graphical  Design  Example 

For  the  purposes  of  illustration  Figure  1  is  used  to  evaluate  the  yield 
pressure,  py,  for  a  sample  motor  case  (Refs  2&4)  in  terms  of  X  and 
6  (see  Tab  3)  derived  from  the  details  given  in  Table  2  .  The  design 
procedure  is  as  follows: 

a  .  Using  the  left  hand  side  of  Figure  1  construct  a  straight 

line  to  pass  through  the  derived  values  of  ft*  -  0 .22  ang 
B  “  1.82.  Mark  the  intersection  of  this  line  with  the  ogw 
scale  . 

b  .  Using  the  lower  right  hand  portion  of  the  nomogram 

construct  a  second  straight  line  to  pass  through  g  -1 .82  and 
X  ■  0  .5  .  Mark  the  intersection  of  this  line  with  the  4 
scale  . 

c.  Project  two  lines,  the  first  horizontally  and  the  second 
vertically,  across  the  ogw  -  4  grid  from  the  intersections 
obtained  in  a  .  and  b  .  respectively  . 

d  .  Evaluate  the  py  value  corresponding  to  the  intersection 

of  the  line  drawn  in  c  •  and  then  py  from  Equation  9  • 

For  the  example  considered.  Figure  1  gives  a  yield  pressure  value  of 
67  .2MPa  .  This  is  in  good  agreement  with  the  numerical  value  of  68  .5MPa  . 

2  .2  Optimisation  Analysis 

2  .2  .1  Theory 

In  the  design  of  rocket  motor  cases  it  is  a  requirement  that  the  structure 
is  manufactured  to  its  optimum  strength  configuration.  For  overwound 
motor  cases  this  configuration  will  clearly  be  a  function  of  the  pre- 
strain  induced  during  winding,  the  motor  rigidity  resulting  from  the 
number  of  applied  fibre  layers  and  the  mechanical  properties  of  the 
constituent  materials  .  A  previous  analysis  (Ref  2)  has  shown  that  there 
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exists  a  maximum  first  yield  pressure  which  is  independent  of  the  winding 
strain,  the  number  of  applied  fibre  layers  and  the  physical  properties  of 
the  fibre  . 

For  the  von  Mises  yield  criterion  this  maximum  is  given  by  Equation  14  . 
The  maximum  is  independent  of  ew  end  n,  (see  Ref  2)  but  the  number  of 
layers  required  to  achieve  it  is  inversely  related  to  the  winding  6train, 
and  vice  versa.  For  a  particular  motor  (Ref  2)  this  trend  is  depicted  by 
the  dotted  line  in  Figure  2  .  For  the  von  Mises  yield  criterion  the 
optimum  n  values  for  a  prescribed  winding  strain  e^,  and  vice  versa,  is 
obtained  (Ref  2)  as  the  solution  to  the  equation 


{ 


4R2 

R2-R  2 
e 


-1+An(2vc-1)} 


_1 _ 

1+An 


/3ewEt;AF  j  (A,n) 


(25) 


For  a  prescribed  n  value  Equation  25  can  be  easily  solved  for  the  optimum 
winding  strain,  but  the  determination  of  the  optimum  number  of  fibre 
layers  for  a  prescribed  winding  strain  would  generally  require  computer 
assistance.  This  problem  can  be  resolved,  however,  if  Equation  25  is 
used  as  the  basis  for  an  optimisation  nomogram. 

As  in  Section  2.1.1,  Equation  25  can  be  simplified  if  a  non-dimensional 
approach  is  adopted.  Using  Equations  7  to  16,  Equation  25  modifies  to: 

K  !>  BS  ■ 

(26) 


where,  as  previously,  vc  »  1/3. 

The  above  equation  can  now  be  written  more  conveniently  as 
fl(fci)  +  f2(X)f3(B)  -  f4(8) 

where 

fl<*W>  -  -/3*w 
f2(X)  -  4/X  -  1 
f 3( 8)  -  l/{( l+8)ln( 1+8)} 
f 4(B)  -  8/ {3( l+8)ln(  1+8)} 


(27) 


(28) 


The  resulting  nomogram  consists  of  two  graduated  straight  loci  representing 
tkr  and  X,  and  a  graduated  8  locus  which  takes  the  form  of  a  curve  .  The 
nomogram  is  depicted  in  Figure  3  with  constructional  details  omitted  . 
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2  .2  .2  Graphical  Design  Example 

For  the  purposes  of  illustration  Figure  3  is  used  to  determine  the  optimum 
winding  strain  for  the  Table  2  motor  in  terms  of  the  derived  variables  8 
and  X  (see  Tab  3).  The  procedure  is  as  follows: 

a  .  Draw  a  straight  line  to  pass  through  the  values  of  1  .82 

on  the  6  scale  and  0  .5  on  the  X  scale  . 

b.  Determine  the  fa,  value  (1.26)  corresponding  to  the 

intersection  of  the  a  .  line  and  the  graduated  ^  locus  . 

Once  ^  has  been  determined,  the  optimum  winding  strain  follows  from  a 
rearrangement  of  Equation  10.  In  this  example  Figure  3  gives  an  optimum 
winding  strain  of  8  .4  x  10-3  .  This  is  in  good  agreement  with  the  numerical 
value  of  8  .5  x  10-3.  The  maximum  yield  pressure  of  92  .OMPa  corresponding 
to  this  optimum  is  calculated  from  Equation  14  .  This  is  36  .9Z  greater 
than  the  yield  pressure  calculated  in  Section  2  .1  .2  . 

In  the  above  analysis  the  optimum  winding  strain  has  been  derived 
theoretically  .  No  attempt  has  been  made  to  assess  whether  the  case  could 
support  the  initial  stress  without  yielding  .  This  can  be  easily  checked 
using  Figure  4  where  the  winding  limit  condition,  Ofa,  “  1  »  ^,ln(l+$) 

(see  Eqn  11),  for  the  variable  ranges  given  in  Table  1  is  plotted.  If  the 
resulting  (\(n,,  8)  coordinate  value  lies  to  the  right  of  the  line  the 
case  has  deformed  plastically  during  winding  and  the  foregoing  analysis 
is  void.  Furthermore,  an  optimised  configuration  cannot  be  achieved,  a 
condition  associated  with  large  throat  diameters  and  small  inertia  loads  . 

3  VISCOEIASTIC  ANALYSIS 

In  the  preceding  section  all  modulus  terms  were  linear  elastic,  and  the 
effects  of  fibre  viscoelastic  relaxation  were  neglected.  As  a  consequence, 
closed-form  analytical  stress  solutions  were  obtained  for  the  fibres  and 
case  from  which  the  required  nomograms  were  constructed. 

For  short  times,  ie  directly  after  winding,  recent  work  has  shown  that 
these  equations  provide  a  simple  and  reliable  method  for  the  elastic 
stress  analysis  of  overwound  rocket  motor  cases  .  For  longer  times, 
however,  these  solutions  may  be  in  error  .  Fibre  viscoelastic  relaxation 
may  have  a  significant  effect  on  the  case  stress  Immediately  before  firing 
and,  as  a  result,  the  first  yield  pressure.  Unfortunately,  for  all  but 
the  most  trivial  of  examples,  numerical  methods  are  required  to  solve  the 
resulting  viscoelastic  equations  (Ref  4;  •  Nevertheless,  analytical 
solutions  have  been  obtained  for  the  fully  relaxed  viscoelastic  condition 
(Ref  4)  .  Since  these  solutions  are  independent  of  the  detailed  nature  of 
the  relaxation  spectrum  (viscoelastic  effects  are  introduced  only  through 
the  long-term  fibre  relaxation  modulus  Ef(*))  the  resulting  equations  are 
valid  for  all  viscoelastic  problems  pertaining  to  this  study. 

Nomograms  for  the  fully  relaxed  condition  are  therefore  developed  in  this 
section.  For  intermediate  times  numerical  methods  (Ref  4)  must  be 
employed  . 
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3  .1  Yield  Initiation  Predictions 

3  .1  .1  Theory 

For  rocket  motor  cases  subjected  to  Instantaneous  pressure  loads  held  for 
Infinitesimally  small  periods  of  time  the  effects  of  fibre  viscoelastic 
relaxation  on  the  pressure-induced  strains  are  usually  neglected.  Hence, 
the  pressurisation  analysis  reduces  to  an  elastic  problem,  with  the 
pressure-induced  stresses  and  strains  related  through  the  conventional 
Hookean  equations  • 

Omitting  details,  and  using  the  notation  previously  defined,  the  long¬ 
term  case  stress  n'(<»)  on  pressurisation  following  complete 

viscoelastic  relaxation  in  the  fibre  can  be  deduced  from  Reference  A  . 


(n) 


J0cp 


(“)  =  -EcewAF1(A,n){ 


A(-)n+Ef(-)/Ef(0) 
1+A(®)n 


+  PR 
t  c ( 1+An ) 


nv„A  R  2 

(l-_§-  )} 


where 


A(«) 


Ef(°,Hf 

Ectc 


(29) 


(30) 


As  before,  it  is  convenient  to  non-dimensionalise  Equation  29  using 
Equations  7  to  16  to  give 


-  (n) 

CT0cp 


(-) 


30w 


(-)  + 


2p 


73X0+8) 


{2+|A} 


(31) 


where 


S0cP(n)<->  “  °0cp(n)<->/oy 

08w(")  ”  *„lnO+B)  rl±£i 

— +6 
a 

a  -  Ef(-)/Ef(0) 


(32) 


(33) 


(34) 

and  ^  is  defined  by  Equation  20.  As  before  vc  ■  1/3.  The  corresponding 
non-dimensional  component  of  axial  stress  is  again  given  by  Equation  18 
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and  since  the  case  remains  linear  elastic  the  yield  criterion  is  given, 
as  before,  by  Equation  19  with  the  quantity  ^n'  appearing  in  that 
equation  replaced  by  Ogcp  n'(")  •  Equations  187P19  and  31  provide  the 


— > - - -  “■>  “ecp  '  ' 

basis  for  nomogram  construction  . 


Equations  18  and  31  can  now  be  used  to  eliminate  oz  and  Often  n^(“) 
respectively  from  the  von  tiises  yield  criterion.  Equation  l?7Pin  favour 
of  p«  in  order  to  predict  the  pressure  for  the  onset  of  yield  .  The 
resulting  equation  is: 


py2  A{i+42-4)  +  Py  73{"2°6w(‘“)*+°ew 
+  09w(“)2  -1-0 


(35) 


On  solution,  the  positive  root  for  py  is  then  substituted  in  Equation  9 
to  obtain  the  yield  pressure  py  . 

The  graphical  representation  of  Equation  35  now  follows  as  for  Equation 
that  is,  the  variation  of  ag*/")  with  4  is  plotted  for  various  py  values 
to  give  the  plot  depicted  in  the  central  portion  of  Figure  5  . 

Since  the  above  analysis  has  been  derived  in  terms  of  the  variable 
combinations  og^**)  and  4,  it  is  necessary,  as  before,  to  construct 
additional  nomograms  so  that  the  required  design  chart  can  be  used  in 
conjunction  with  the  fundamental  variables  B,  o  and  X.  Although 
the  graphical  determination  of  4  readily  follows  from  the  Equation  23 
nomogram,  the  graphical  determination  of  the  effects  of  complete  fibre 
relaxation  on  the  case  pre-stress,  08W<“),  cannot  be  established  from 
the  ’z'  type  nomogram  shown  to  the  left  of  the  central  plot  in  Figure 
1,  since,  when  viscoelastic  effects  are  introduced,  the  case  pre-stress 
becomes  a  function  of  both  the  long-  and  short-term  fibre  modulus  values 
(see  Eqn  33)  .  ^Under  these  conditions,  the  effects  of  complete  fibre 
relaxation  on  oewC9)  may  be  graphically  determined  from  a  ’proportional' 
type  nomogram  (Ref  6)  of  the  form 


fl(o,8)  _  f3(4„) 

f2(B)  "  f^OewC-)) 

where  the  functions  represent 
fl(a.B)  "  1/a  +  «  -  y 
(2(B)  -  (1+B)ln(l+B) 
f3(*w)  -  *w 
/4(»9w(*))  "  aew(*) 


(36) 


(37) 
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As  before,  the  three  Independent  nomograms  can  be  combined  through  their 
common  loci  •'  The  final  design  chart  is  shown  in  Figure  5  with  constructional 
details  omitted. 

In  contrast  to  Figure  1  the  design  process  embodied  In  Figure  5  Is  now 
In  terms  of  five  fundamental  variables  tw>  X,  py,  a  and  B*  However,  the 
variable  6  Is  again  represented  by  two  of  the  graduated  loci  in  the  design 
chart.  (Furthermore,  y  is  a  function  of  BO  It  follows  therefore,  as 
before,  that  6  cannot  be  'read'  directly  from  the  design  chart  in  terms  of 
the  variables  py,  X  and  a  and  that  a  trial-and-error  solution  procedure 
is  again  required,  ffevertheless ,  this  problem  can  be  resolved  if  the  fully 
relaxed  case  pre-stress,  oqwC*"),  is  specified  together  with  X,  py  and  a  so 
that  the  solution  procedure  follows  in  a  manner  similar  to  that  described 
in  Section  2  .1  .1  . 

3  .1  .2  Graphical  Design  Example 

To  illustrate  the  use  of  Figure  5  the  fully  relaxed  viscoelastic  yield 
pressure  value  is  calculated  for  the  Table  2  motor  using  the  derived 
variables  given  in  Table  3.  However,  as  Figures  1  and  5  are  Identical, 
except  for  the  case  pre-stress  nomogram,  only  the  use  of  the  case  pre¬ 
stress  nomogram  will  be  discussed  here.  The  procedure  is  as  follows: 

a  .  Using  the  left-hand  side  of  the  nomogram,  construct  an 

index  line  to  pass  through  1  .82  on  the  B  scale  and  3  .82  on  the 
y  scale  .  Mark  the  intersection  of  this  line  with  the  ungraduated 
diagonal  . 

b  .  From  the  intersection  established  in  a  .  draw  a  second 

straight  line  to  cut  the  Bw  scale  at  0.22.  Extend  to  cut  the 
ogu(“)  scale  at  0.17.  This  is  the  non-dimensional  fully 
relaxed  viscoelastic  case  pre-stress  . 

c.  Having  established  use  steps  b  .  to  d  .  of  Section 

2  .1 .2  to  calculate  the  yield  pressure  as  62  .98MPa  •  This  compares 
favourably  with  the  numerical  value  of  62  .4MPa  . 

The  effects  of  viscoelastic  relaxation  can  now  be  assessed.  By  comparing 
the  yield  pressure  values  obtained  from  the  elastic  and  fully  relaxed 
viscoelastic  conditions  (Figs  165)  for  the  Table  2  motor  case,  ie  67 .2MPa 
and  62  .9MPa  respectively,  it  can  be  seen  that  the  yield  pressure  has 
decreased  by  approximately  6  .5X .  It  should  be  noted ,  however ,  that  the 
yield  pressure  may  not  always  decrease  .  As  has  been  reported  elsewhere 
(Ref  4),  depending  on  the  initial  winding  configuration,  the  yield  pressure 
may  Increase,  decrease  or  Increase  to  a  maximum  and  then  decrease  with  fibre 
relaxation.  Nevertheless,  the  maximum  attainable  yield  pressure  will 
always  be  given  by  Equation  14  . 

An  insight  into  which  of  these  three  yield  pressure  trends  a  motor  will 
exhibit  can  be  obtained  from  a  consideration  of  the  ogw  -  6  grids  in 
Figures  1  and  5  •  If  the  relevant  point  lies  to  the  right  of  the  ^  *  1 
curve  in  Figure  1  the  yield  pressure  will  always  decrease  with  fibre 
relaxation,  whereas  the  opposite  is  true  if  the  intersection  lies  to  the 
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left  of  the  p„  »  1  curve  In  Figure  S.  But,  If  the  Intersection  lies  to 
the  left  of  tne  py  ■  1  curve  In  Figure  1  and  then  lies  to  the  right  of 
the  py  •  1  curve  in  Figure  5,  the  yield  pressure  will  have  increased  to  the 
maximum  and  then  decreased  with  fibre  relaxation.  For  design  purposes 
the  latter  trend  Is  preferable,  provided  the  elastic  and  fully  relaxed 
viscoelastic  yield  pressure  values  are  equal.  Otherwise,  depending  upon 
Initial  winding  conditions,  the  yield  pressure  value  at  either  the  elastic 
or  the  fully  relaxed  viscoelastic  condition  may  take  an  unacceptably  low 
value  . 

3  .2  Optimisation  Analysis 

3  .2  .1  Theory 

In  the  preceding  section  it  was  shown  that  the  yield  pressure  depends  on 
fibre  relaxation  .  It  follows  therefore  that  the  elastic  optimisation 
nomogram  (Fig  3)  is  not  valid  when  viscoelastic  relaxation  occurs  . 

This  is  illustrated  in  Figure  2,  where  the  first  yield  pressure  values 
are  plotted  against  the  number  of  applied  fibre  layers  for  a  particular 
motor  (Ref  4)  for  the  elastic  and  fully  relaxed  conditions  .  It  can  be 
seen  that  fibre  relaxation  is  accompanied  by  a  shifting  of  the  yield 
pressure  curve  so  that  the  number  of  fibre  layers  required  for  the 
fully  relaxed  optimum  configuration  is  greater  than  that  required  for  the 
elastic  analysis  .  Furthermore,  for  the  elastic  optimised  motor  the  yield 
pressure  in  the  fully  relaxed  condition  will  be  lower  than  that  given  by 
Equation  14. 

From  a  previous  viscoelastic  analysis  (Ref  4)  it  can  be  shown,  after  some 
manipulation,  that  when  viscoelastic  effects  are  introduced  the  optimisation 
equation  for  the  fully  relaxed  condition  becomes 


4R2 

R*-Re* 


-1+An(2vc-1)} 


1+An 


✓3twEcAFl(A.n)  {  1+An 


Ef(0) 

Ef(-) 


+An 


(38) 


For  a  prescribed  number  of  fibre  layers  Equation  38  can  be  solved  for  the 
fully  relaxed  optimum  winding  strain  .  Conversely,  the  determination  of 
the  optimum  number  of  fibre  layers,  for  a  prescribed  winding  strain, 
would  generally  require  computer  assistance .  It  is  therefore  necessary 
to  represent  Equation  38  graphically . 

Using  Equations  7  to  16  and  Equations  32  to  34  and  setting  vc  _  1/3, 
Equation  38  reduces  to  the  following  non-dimensional  form: 


d-1-  ») -jlj  .  W),^, 


(39) 
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which  on  subsequent  rearrangement  corresponds  to  an  alignment  nomogran 
(Ref  6)  of  the  form 


fj(a)  +  f2(*w)*3(X.8)  “  *4(8) 

where  the  functions  represent 
fl(a)  ■  l/o 
*2(*w)  “  -*w 

f,(X,B)  -  {/3 ( 1+8) 2ln(  1+8 )}/{£- 1-8} 
3  X  3 

*4(6)  -  -8 


(40) 


(41) 


This  is  analogous  to  Equation  27  .  The  resulting  nomogram  will  therefore 
consist  of  two  graduated  parallel  loci  for  and  a  and  a  series  of 
graduated  curved  8  loci,  each  for  a  prescribed  X  value. 

In  principle  the  resulting  nomogram  will  graphically  represent  Equation 
39.  In  practice,  however,  experience  has  shown  that  unacceptable  reading 
errors  can  occur  for  extreme  variable  combinations  .  To  resolve  this 
problem  it  is  necessary  to  transform  the  Equation  40  nomogram  into  an 
alternative  form  so  that  the  resulting  design  chart  possesses  greater 
practical  utility. 


Using  determinant  theory  (Ref  7)  and  omitting  the  details  it  can  be  shown 
(see  Annex  A)  that  the  required  nomogram  can  be  obtained  from  the  basic 
determinant 


-l/*l(o)  0  1 

0  -l/f2(<M  1 

-1/f 4( 8)  -f3(X,B)/f4(8)  1 


-  0 


(42) 


This  corresponds  to  a  nomogram  which  consists  of  a  locus  for  a 
coincident  with  the  x-axis,  a  locus  of  tw  coincident  with  the  y-axis 
and  a  series  of  8  curves,  each  curve  representing  a  unique  X  value. 
Furthermore,  the  m  and  ♦w  axes  are  orthogonal .  The  resulting  nomogram  is 
presented  in  Figure  6  with  the  constructional  details  omitted* 

Alternative  transformations  of  Equation  40  can  be  derived  leading  to 
different  determinant  forms,  but  these  have  not  led  to  significantly 
better  nomographic  representations. 
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3  -2  .2  Graphical  Design  Example 

For  the  purposes  of  Illustration,  Figure  6  Is  used  to  determine  the 
winding  strain  necessary  to  obtain  the  optimised  configuration  at  the 
fully  relaxed  condition  for  the  Table  2  motor  In  terms  of  the  derived 
variables  given  in  Table  3.  The  procedure  is  as  follows: 

a.  Locate  a  point  within  the  X  -  6  grid  corresponding  to  the 

intersection  of  8  “  1  -82  and  X  “  0  .5  . 

b  .  From  this  point  draw  a  straight  line  to  pass  through 

0.5  on  the  a  scale  to  obtain  the  required  ^  value. 

For  the  example  considered.  Figure  6  gives  a  winding  strain  of  1  .17  x 
10~2,  as  compared  to  the  numerical  value  of  1.15  x  10“2.  These  values  are 
in  good  agreement  .  If  the  elastic  and  fully  relaxed  viscoelastic  optimum 
winding  strains  are  compared  (8.46xlO-3  and  1  .17xl0-2  respectively)  it  can 
be  seen  that  the  latter  is  38  .32  greater  than  the  former  .  Consequently, 
it  may  not  be  possible  to  obtain  a  fully  relaxed  optimised  configuration; 
the  pre-stress  induced  during  winding  may  exceed  the  elastic  limits  of 
the  case  material.  This  condition  can,  as  previously,  be  easily  checked 
using  Figure  4  .  If  the  relevant  point  lies  to  the  right  of  the  winding 
limit  condition  line,  ie  a^,  »  1  •  fwlnO+B)  (see  Etln  H)»  tlle  case 
has  deformed  plastically  during  winding  and  the  foregoing  analysis  is 
void  . 

4  MAXIMUM  FIBRE  STRESS 

4 .1  Theory 

In  the  foregoing  sections,  graphical  design  procedures  were  developed  to 
determine  the  maximum  permissible  applied  pressure  and  the  optimised 
winding  configuration  for  both  the  initial  elastic  and  fully  relaxed 
conditions  .  In  these  analyses  the  staxlmum  pressure  calculation  has  been 
based  on  the  pressure  for  first  yield  in  the  case  material,  and  the 
magnitude  of  the  fibre  stresses  associated  with  this  pressure  haB  not 
been  studied  •  In  the  present  application  it  is  likely  that  this  approach 
will  provide  a  satisfactory  design  solution,  since  the  ultimate  tensile 
strength  of  the  fibre  overwind  is  very  much  greater  than  the  yield  stress 
of  the  case,  fevertheless,  care  must  be  taken  at  the  exploratory  stages 
of  the  design  to  ensure  that  the  fibre  stresses  do  not  exceed  prescribed 
safety  values  . 

A  previous  study  (Ref  2)  has  shown  that  the  maximum  fibre  stress  occurs 
in  the  outermost  fibre  layer  directly  after  winding.  It  follows  therefore 
from  Equation  1  that  the  maximum  fibre  stress  09 fmax  18  given  approxi¬ 
mately  by 


o 


9f  max” 


*wEf{l- 


1+An 


+  P** 

tf(l+An) 


n-^i- 


(43) 
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where,  as  previously,  vc  •  1/3.  It  should  be  noted  that  the  above 
equation  is  also  valid  for  the  viscoelastic  fibre  with  Ej  replaced  by 
Ef(0). 

Again  it  is  convenient  to  adopt  a  non-dimensional  approach.  Using 
Equations  7  to  16  it  can  be  shown  that  Equation  43  reduces  to: 


'w  «'V  If  - 

<«> 

where  5  denotes  the  modulus  ratio  Ef/Ec  and  Ogf^x  is  the  non- 
dimensional  maximum  fibre  stress  Ogfj^x/oy  .  In  Equation  44 
the  approximation 


1  -  A  ♦  1 
It- An 


(45) 


has  been  used . 

One  of  the  most  simple  techniques  (Ref  7)  for  dealing  with  equations  of  the 
above  form  is  to  Introduce  an  additional  auxiliary  variable  as  follows  . 

°0f  max  “  S<*w+,c> 

(46) 

where 


K 


2  P 
75  TwST 


(2  -  1) 
a  T 


(47) 


Both  of  these  equations  can  now  be  readily  represented  by  grid  nomograms. 
Equation  46  by  an  alignment  (Ref  6)  nomogram  and  Equation  47  by  a  'r'  type 
(Ref  6)  nomogram.  The  complete  nomogram  is  then  obtained  by  combining  these 
two  Independent  nomograms  in  such  a  way  that  the  locus  k  is  common  to 
both.  It  should  be  noted,  however,  that  if  a  conventional  construction 
technique  is  used  for  the  Equation  46  nomogram,  the  desired  grid  will  not 
be  obtained,  since  the  effect  of  varying  A  is  merely  to  adjust  the 
gradations  on  the  6  scale,  not  to  alter  its  spatial  position  relative  to 
the  nomogram  as  a  whole.  Nevertheless,  a  pseudo-grid  can  be  constructed 
whereby  the  6  scales,  for  each  unique  A  value,  are  projected  horizontally 
onto  a  pseudo-scale  set  off  to  the  right-hand  side  of  the  actual  0  locus  . 

By  passing  loci  through  common  variable  values  on  these  scales  the  required 
grid  is  formed  .  The  ordinate  of  any  point  in  the  grid  is  therefore  projected 
horizontally  across  to  the  locus  which  forms  the  left  hand  edge  of  the  6-A 
grid.  This  process  can,  for  a  known  intersection  of  the  left-hand  edge 
of  the  0-A  grid,  be  reversed  if  necessary  to  obtain  A  for  a  known  0, 
and  vice  versa . 
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The  complete  nomogram  Is  depicted  In  Figure  7  with  constructional  details 
omitted  .  The  common  locus  ic  Is  represented  by  the  ungraduated  locus 
located  at  the  centre  of  the  nomogram  • 

4  .2  Graphical  Design  Example 

For  the  purposes  of  illustration  the  maximum  fibre  stress  is  calculated 
for  the  Table  2  motor  using  the  derived  variables  given  in  Table  3  • 

The  applied  internal  pressure  is  taken  as  that  for  first  yield  in  the 
case  material  assuming  elastic  conditions  (see  Section  2  .1  .2)  .  The 
procedure  is  as  follows: 

a.  Using  the  right-hand  side  of  the  nomogram  determine  the 

intersection  of  the  X  ■  0 .5  and  8  “  1-82,  Interpolating  as  necessary. 
Project  the  ordinate  of  this  point  horizontally  to  intersect  the 
left-hand  edge  of  the  B-X  grid. 

b  .  Construct  an  index  line  to  pass  through  the  final 

Intersection  obtained  in  a.  and  the  value  0.73  on  the  p  scale  to 
cut  the  locus  positioned  at  the  centre  of  the  nomogram. 

c.  Construct  a  second  index  line,  this  time  to  pass  through 

the  grid  point  corresponding  to  5  “  1.77  and  •  0.22  and  the 
intersection  of  the  k  locus  determined  in  b  .  to  obtain  the  required 
°efmax  • 

The  maximum  fibre  stress  of  llOOMPa  then  follows  from  the  simple  relationship 
°6f max  “  ^flfmax  x  °y  •  Thi®  in  good  agreement  with  the  numerical 
value  of  1069MPa  and  consistent  with  the  approximations  made  in  Equation 
44.  The  UTS  of  the  fibre  used  in  this  example  (Tab  2)  is  2400KPa  . 

5  CONCLUSIONS 

Charts  have  been  presented  for  use  in  the  design  of  a  metallic  cylinder 
reinforced  with  a  viscoelastic  fibre  overwind.  Topics  considered  Include: 

a.  Determination  of  the  pressure  for  first  yield  in  the 

metallic  case  at  short-  and  long-times  . 

b  .  Determination  of  the  optimised  configuration  for  short- 

and  long-times  • 

c  .  Determination  of  the  maximum  fibre  stress  . 

Examples  have  been  given  with  each  nomogram  for  a  particular  component; 
derived  results  compare  favourably  with  exact  numerical  results. 
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7  NOMENCLATURE 

Variables 

A,A(«)  defined  by  Equations  3  and  30  respectively 
E,E(«)  short-  and  long-term  Young's  modulus  values  respectively 
Fi(A,n)  defined  by  Equation  4 
n  number  of  fibre  layers 

p  Internal  pressure 

R.Rg  internal  case  radius  and  effective  nozzle  radius  respectively 
t  thickness 

e  strain 

v  Poisson's  ratio 

a  stress 

Derived  Variables 

p  applied  /maximum  pressure  ratio  (see  Eqn  9) 

a  fibre  viscoelastic  modulus  parameter  (see  Eqn  34) 

6  stiffness  parameter  (see  Eqn  7) 

y  stiffness-fibre  modulus  parameter  (see  Eqn  37) 

6  constituent  materials  modulus  ratio  Ef/Ec 

*  auxiliary  variable  (see  Eqn  47) 

X  radius  parameter  (see  Eqn  8) 

4  auxiliary  variable  (see  Eqn  20) 

♦w  non-dimensional  winding  parameter  (see  Eqn  10) 

Subscript 

f,c  fibre  and  case  respectively 

max  maximum 

opt  optimum 

p  due  to  pressurlsation 

uts  ultimate  tensile  strength 

w  winding 

y  at  first  yield 

8,z  circumferential  and  axial  directions  respectively 

Superscript 

(i,n)  i'th  fibre  of  n  fibre  layers 

(n)  n  fibre  layers 

Accent 

a  non-dimensional 
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TABLE  1  Topical  Ranges  of  fon-Dimensional  Variables 


Values 
Min  May 


tbn-dimenslonal  winding  parameter 

Jw 

0 

1  .5 

”  case  pre-stress 

Offw 

0 

1  .0 

Applied/maximum  pressure  ratio 

A  “ 

P 

0 

1  .0 

Radius  parameter 

X 

0 

1 .0 

Stiffness  parameter 

8 

0 

5.0 

Fibre  viscoelastic  modulus  ratio 

a 

0 

1  .0 

Stiffness-fibre  modulus  parameter 

y 

0 

15  .0 

Constituent  materials  modulus  ratio 

s 

0 

4  .0 

Ifon-dimensional  maximum  fibre  stress 

°0fmax 

0 

15  .0 

TABLE  2  Sample  Motor  Case  Details 


Fibre  moduli 

Short  term 

Ef(0) 

124GPa 

Long  term 

Ef(-) 

62GPa 

Fibre  UTS 

Outs 

2400MPa 

Case  Young's  modulus 

Ec 

70GPa 

Case  yield  stress 

Oy 

470MPa 

Thicknesses 

fibre 

tf 

0  .1mm 

case 

tc 

1  .95mm 

Radii 

Motor 

R 

46  .Omra 

Effective 

32  .53mm 

Winding  strain 

*w 

1  .5  x  10“3 

Mimber  of  layers 

n 

20 

TABLE  3  Derived  Variables  for  Table 

Iton-dlmenslonal  winding  parameter 

2  Motor 

0.22 

Radius  parameter 

X 

0.5 

Stiffness  parameter 

B 

1 .82 

Fibre  viscoelastic  modulus  ratio 

a 

0.5 

Stiffness-fibre  modulus  parameter 

Y 

3  4(2 

Constituent  materials  ratio 

6 

1.77 
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APPLIED/MAXIMUM  PRESSURE  RATIO 
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LAYERS  OF  REINFORCEMENT 
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FIG.  3 
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FIG.  3  NOMOGRAM  RELATING  NON-DIMENSIONAL  DESIGN  PARAMETERS 
FOR  MAXIMUM  FIRST-YIELD  PRESSURE:  ELASTIC  OVERWIND 


UNLIMITED 


FIG.  6 


a  FIBRE  VISCOELASTIC  MODULUS  PARAMETER  B  STIFFNESS  PARAMETER 


FIG.  6  NOMOGRAM  RELATING  NON-DIMENSIONAL  DESIGN 
PARAMETERS  FOR  MAXIMUM  FIRST-YIELD  PRESSURE: 
FULLY  RELAXED  VISCOELASTIC  OVERWIND 
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FIG.  7  NOMOGRAM  RELATING  NON-DIMENSIONAL  MAXIMUM  FIBRE 
STRESS  TO  NON-DIMENSIONAL  DESIGN  PARAMETERS 
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AN1CX  A 


DETERMINANT  THEORY  FOR  NOMOGRAM  COKSTRUCTION 

If  the  coordinate  point  x^.y^  Is  represented  parametrically  In  terms  of 
a,  le 


*1  “  xi(a)  yi  "  yi(a) 


(Al) 


a  curve  is  defined  in  the  x-y  plane  .  Furthermore,  if  a  is  specified  the 
corresponding  point  on  the  curve  is  determined  .  Consequently  the  curve 
can  be  graduated  in  terms  of  the  parameter  a  • 

Similarly,  two  further  curves  can  be  defined  parametrically  by  the 
equations 


X2  “  *2(b) 

y2  ■  y2<b> 

X3  -  X3(c) 

y3  ■  y3(c> 

(A2) 

(A3) 


and  graduated  in  terms  of  the  parameters  b  and  c  respectively  • 

If  the  three  coordinate  points  x^.y^,  X2»y2  and  X3,y3  are  to  be  co- 
linear  it  is  a  requirement  that  the  gradient  of  the  line  joining  x^.y^ 
and  X2,y2  la  the  same  as  that  of  the  line  joining  x^.y^  and  X3,y3  le 


yz~y\  .  ?3-ri 
x2-xi  x3-xi 

(A4) 

which,  on  expanding,  and  using  the  parametric  relationships  defined  above 
gives 

Xl(a){y2(t>)-y3(0>  +  *2(bHy3(c)-yi<a>>  +  *3(0  {yi(a)-y2(b)  )  “  0 

(A5) 


This  equation  can  also  be  expressed  in  determinant  form  (Ref  7)  as 


xj  (a) 

yi(0 

x20>) 

y2(b) 

X3(c) 

y3(c> 

(A6) 


The  characteristics  of  the  determinant  in  Equation  A6  are: 

a.  Its  value  is  aero, 

b  •'  each  row  contains  one  parameter  only, 

c  •  each  element  in  the  last  row  is  equal  to  unity  • 


Al 
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It  is  now  necessary  to  regard  the  parameters  a,  b  and  c  in  Equation  A6 
as  the  variables  of  a  three-variable  equation.  If  (and  only  if)  such  an 
equation  can  be  reduced  to  a  determinant  of  the  form  of  Equation  A6,  it 
follows  that  a  solution  to  the  equation  can  be  obtained  from  the  points 
of  intersection  of  any  straight  line  with  the  three  graduated  curves 
given  by  Equations  Al,  A2  and  A3,  ie  a  nomogram  can  be  constructed.  The 
form  of  the  curves  in  the  nomogram  is  determined  by  the  elements  of  the 
corresponding  row  in  the  determinant  (Ref  7)  . 

If,  in  Equation  40,  f(a)  is  represented  by  x  and  f2('kj)  by  y»  iben  the 
following  three  simultaneous  equations  may  be  written 


x  -  fi(o) 

y  “ 

x  +  f  3(  X,  8)y  -  f  4(  8) 
ie 

x  +  Oy  -  fi(a)  “  0 
Ox  +  y  -  f2(*w>  “  0 
x  +  f  3(  X,  B)y  -  f  4(  0)  “  0 

The  condition  for  a  non-trlvial  solution  of  Equation  A8  is 


(A7) 


(A8) 


0 

1 

f 


-f  l(o) 

-f4<B) 


(A9) 


This  determinant  is  not  of  the  same  form  as  that  of  Equation  A6,  but  by 
applying  standard  determinant  operations  (Ref  7),  the  following  transformed 
version  can  be  derived 


-1/f  l(a) 
0 

-l/f4(B) 


0  1 

-l/f2<*w>  1 

-f  3(X,8)/f4(6)  1 


(A10) 


where  X  is  considered  to  be  a  constant  •  Equations  A6  and  A10  are  now  of 
the  same  form.  Consequently,  it  can  be  concluded  that  Equation  40  can  be 
represented  nomographlcally  as  detailed  in  the  main  text  following  Equation 


42. 


Further  transformations  of  Equation  A10  can  also  be  performed  leading  to 
different  determinant  forms,  but  these  have  not  led  to  significantly  better 
nomographic  representation  of  Equation  40 . 
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